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Conjectured Exact Locations of Dynamical Transition Points for 
the ±J Ising Spin Glass Model 

Chiaki Yamaguchi 
Kosugichou 1-359, Kawasaki 211-0063, Japan 

The conjectured exact locations of the dynamical transition points for the ± J Ising spin 
glass model are theoretically shown based on a percolation theory and conjectures. The 
dynamical transition is a transition for the time evolution of the distance between two spin 
configurations. The distance is called the damage or the Hamming distance. The conjectured 
exact locations of the dynamical transition points are obtained by using the values of the 
threshold fractions of the random bond percolation problem. The present results are obtained 
as locations of points on the Nishimori line which is a special line on the phase diagram. 
We obtain Td = 2/\n{z/z — 2) and p D = z/2(z — 1) for the Bethe lattice, Td oo and 
p D — > 1/2 for the infinite-range model, Td = 2/ In 3 and p D — 3/4 for the square lattice, 
T D ~ 3.9347 and p D ~ 0.62441 for the simple cubic lattice, T D ~ 6.191 and p D ~ 0.5801 
for the 4-dimensional hypercubic lattice, and Td = 2/ln[l + 2 sin(7r/18)/l — 2sin(7r/18)] 
and p D = [1 + 2 sin(7r/18)]/2 for the triangular lattice, when J/ttB = 1, where z is the 
coordination number, J is the strength of the exchange interaction between spins, ks is the 
Boltzmann constant, Td is the temperature at the dynamical transition point, and p D is the 
probability, that the interaction is ferromagnetic, at the dynamical transition point. 

§1. Introduction 

To establish reliable analytical theories of spin glasses has been one of the most 
challenging problems in statistical physics for years .OJSJSJllJDJZHni Our main in- 
terest in this article does not lie directly in the issue of the properties of the phases 
in spin glasses. We instead will concentrate ourselves on the precise determination 
of the structure of phase diagram for dynamical properties of spin glasses. This 
problem is of practical importance for numerical studies, since exact locations of 
transition points greatly facilitate reliable estimates of physical properties around 
the transition points. 

The dynamical transition mentioned in this article is a transition® '^'D®'^ for 
the time evolution of the distance between two spin configurations. The distance 
is called the damage or the Hamming distance. In the distance measurement, the 
time evolution of two spin configurations is studied under the same thermal noise, 
where, as for the initial condition, the two configurations are generally set to have a 
relationship with each other. 

The ± J model is known as one of the Ising spin glass models.^'®'®'®'®'® 1 A line 
on the phase diagram for the ± J model is called the Nishimori line.® Interestingly, 
the internal energy, the upper bound of the specific heat and so forth are exactly 
calculated on the line.®'®'®'®'®'® In addition, the internal energy does not depend 
on any lattice shape and instead depends on the number of nearest-neighbor pairs 
in the whole system. The location of the multicritical point is conjectured on the 
Nishimori line of the square lattice, and the conjectured location is in good agreement 
with the results of other numerical estimates.^ The present results are also obtained 
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as locations of points on the Nishimori line. 

We use two conjectures in this article. One of the conjectures has already been 
proposed and discussed in Refs. EI))©) and another will be proposed in this article. 

In order to obtain the conjectured exact locations of the dynamical transition 
points, the values of the threshold fractions of the random bond percolation problem 
are used. If a threshold fraction of the random bond percolation problem is calcu- 
lated, one is able to calculate the conjectured location of the dynamical transition 
point in the ±J model by using the present theory. Generally, calculation of the 
threshold fractions of the random bond percolation problem is easier than that of 
the dynamical transition points in the ± J model. Therefore, the present theory can 
be promising in this respect. 

This article is organized as follows. In £J21 the ± J model is explained. In £}3l our 
conjectures are described, and conjectured exact equations are shown. The present 
results by using the obtained equations are given in §H In £j5l the concluding remarks 
of this article are described. 

§2. Model 

The Hamiltonian for the ±J model, H, is given by^'^ 

7~L = — JijSiSj , 

where denotes nearest-neighbor pairs, Si is a state of the spin at the site i, and 
S{ = ±1. Ji j is the strength of the exchange interaction between the spins at the 
sites i and j. The value of Jy is given with a distribution P(Jij). The distribution 
P(Jij) is given by 

P(Aj) /"U ,./ + (1 " P) S JiJt -j , (2-2) 

where J > 0, and 5 is the Kronecker delta, p is the probability that the interaction is 
ferromagnetic, and 1 — p is the probability that the interaction is antiferromagnetic. 

The dynamical transition mentioned in this article is characterized by a dis- 
tance between two spin configurations {Si} and {Si}. It has been pointed out that, 
for the cubic lattice, there are three phases,^ i.e., a high-temperature phase, an 
intermediate phase and a low-temperature phase. In the high-temperature phase, 
the two configurations become identical quickly, so that the distance between them 
vanishes. In the intermediate phase and the low-temperature phase, the distance 
between the two configurations remains a finite in the long-time limit if the system 
size is large enough. In the intermediate phase, the distance between the two config- 
urations does not depend on the initial conditions of the two configurations. In the 
low-temperature phase, the distance between the two configurations depends on the 
initial conditions of the two configurations. We concentrate ourselves on the dynam- 
ical transition between the high-temperature phase and the intermediate phase, and 
especially concentrate ourselves on the location of the dynamical transition point. 

We explain the measurement of the distanceP^ 1 '^ The distance D(i) is given 
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by 

1 N 

D ® = l 5 *(*) " S M)mc}r , (2-3) 

i 

where t is the time, N is the number of sites, ( )mc denotes the sample average 
by the Monte Carlo method, and []r denotes the random configuration average for 
the exchange interactions. As for the initial condition of {Si} and {Si}, {S'j(O)} 
is set to random, and {Sj(0)} = — {S'j(O)} is set, for example. As for the Monte 
Carlo method, a heat-bath method is used, and is as follows. By using the uniform 
pseudo-random number r(t) (0 < r(t) < 1), update of {Si(t)} and {Si(t)} is given 

by 

S ,( t+ ^ ) = ^ IT - F ^L__ _ rW |, (2 . 4) 
^ + ^ ) ^^ __^__- - r(t) }, ( , 5) 



where At is a time step per spin, the summations of the right-hand sides of Eqs. (|2^ 
and (|2-5p are over the nearest-neighbor sites of the site i, (3 = l/kpT, T is the 
temperature, and kp is the Boltzmann constant. Note that, in order to measure 
D(t), the two spin configurations {Si(t)} and {Si(t)} are updated with the same 
pseudo-random number sequence {r(t)}. 

We apply a percolation theory. We use the Fortuin-Kasteleyn (FK) cluster EJ'tHJ 
The FK clusters consist of the FK bonds which are probabilistically put between 
spins. The number of the FK bonds is rigorously related to the internal energyP 1 
We define the probability for putting the FK bond as -Pfk- -Pfk is given 

Pfk = 1 - e -0 J « SiS i-fl J «l . (2-6) 

For calculating [{Pfk)t]r, a gauge transformation is used, where ()t denotes the 
thermal average. The gauge transformation is performed by^ 1 '^ 

Ji,j — > Jij&i&j , Si — > SiCTi , (2-7) 

where Cj = ±1. The gauge transformation has no effect on thermodynamic quanti- 
ties.^ By performing the gauge transformation, the T-L part becomes % — > % and 
the Pfk part becomes Pfk —> Pfk- By using Eq. (|2-2p . the distribution P(Jjj) is 
rewritten aP 

P( ^ )= 2cosh(/ ? pJ) - J "i = ±J - (2 - 8 » 

where j3p is given by 

fr^ln^. (2-9) 
By performing the gauge transformation, the distribution P(Jij) part becomes 



[2cosh(/3 P J)]^ 
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^ (2-10) 

where Np is the number of nearest-neighbor pairs in the whole system. When the 
value of j3p is consistent with the value of the inverse temperature f3, the line for 
P = (3p on the phase diagram is called the Nishimori lineP By using the gauge 
transformation, [(Pfk)t]r on the Nishimori line is obtained as^® 

[(Pfk)t]r = tanh(/3 P J) , (2-11) 

where f3p = 1/kpTp, Tp is the temperature on the Nishimori line. The value of 
[{Pfk)t]r on the Nishimori line does not depend on any lattice shape. 

§3. Conjectures 

The present theory uses two conjectures. If the two conjectures are correct, the 
present theory gives the exact values. Here, we describe the two conjectures and 
show conjectured exact equations derived by using the two conjectures. 

(1) We conjecture 

T FK = T D (3-1) 
for arbitrary lattices, where Tfk is the percolation transition temperature of the FK 
cluster, and Tp> is the dynamical transition temperature for the distance between two 
spin configurations. Two-spin correlations of nearest neighbor pairs are correctly 
estimated by using the FK bonds, while, in the case of the ±J model, two-spin 
correlations of no nearest neighbor pairs are not correctly estimated from simple 
connection of the FK bonds, however, the transition point for the time evolution of 
the distance between two spin configurations may be consistent with the percolation 
transition point of the clusters for the net two-spin correlations by the FK bonds. 
This conjecture for Tfk = Tp, has already been proposed and discussed and the 
confirmations of this conjecture have been numerically performed.®® The previous 
studies indicate that Tfk may exactly agree with T^P'® 

(2) We conjecture 

[{Pfk)t]r = Pc (3-2) 
at the percolation transition point of the FK cluster on the Nishimori line of arbitrary 
lattices, where Pfk is the probability for putting the FK bond between spins, and 
Pc is the threshold fraction of the random bond percolation problem. In the random 
bond percolation problem, bonds for generating clusters are randomly put on the 
edges of the lattice, and one of the clusters is percolated at the threshold fraction 
P c m From Eq. OTTD . the value of [{Pfk)t\ K on the Nishimori line does not 
depend on any lattice shape. From the fact, Eq. (j3-2|) is conjectured. We propose 
this conjecture in this article as a conjecture which may be exact. 
By using Eqs. (f2~TTjh ([3T]) and ([32]), we obtain 

Td = k B Hl+ 2 p c /l-Pc) ' (3 ' 3) 
P D = \{1 + P C ). (3-4) 
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Eqs. ()3-3p and (|3-4p are conjectured exact equations. By using Eqs. (j3-3|) . (|3-4p and 
the value of the threshold fraction Pc of the random bond percolation problem, 
the values of the dynamical transition temperature Td and the dynamical transi- 
tion probability p D are calculated as the location of a point on the Nishimori line. 
The obtained values are conjectured exact values. Note that the dynamical tran- 
sition probability p D is the probability that the interaction is ferromagnetic at the 
dynamical transition point. 

Campbell and Bernardi have derived an equation for the energy E in the ±J 
model and the threshold fraction Pc of the random bond percolation problem, i.e., 
E = JNb{1 — 2Pc/[l — exp(— 2/3J)]} on the assumption of a random active-bond 
spatial distribution.® By applying the energy and the temperature on the Nishimori 
line to this equation, the same equations (Eqs. (|3-3p and (|3-4p ) are obtained]-® where 
the energy on the Nishimori line is — Nb Jtanh(/3p J)P 

In Ref. ITS]) , the exact values of the percolation threshold of the FK cluster are 
also conjectured. The conjecture proposed in this article is more general, so that the 
present study includes finite-dimensional cases that are not mentioned in Refll5p. 

§4. Results 

We show the present results by applying the conjectured exact equations ob- 
tained in £}3l i.e., Eqs. (|3-3p and (|3-4p . to several lattices. The present results are 
obtained as locations of points on the Nishimori line. 

For the Bethe lattice, the threshold fraction Pc of the random bond percolation 
problem is obtained as Pc = l/(z — l)p^ where z is the coordination number. By 
using Eqs. (|3-3p and (|3-4p . we obtain 

- 2J _ z 

D ~ k B Hz/z - 2) ' Pd ~ 20-1) ' 

This is the result for the Bethe lattice. Tjj agrees with the ferromagnetic transition 
temperature for the pure system, TcP3 When z = N — 1 and J — > J/y/N, this 
model becomes the infinite-range model. Then, in the thermodynamic limit, we 
obtain Tjj — > oo and p D — > 1/2 when J/ks = 1. This result for the infinite-range 
model agrees with the previous results in Refs. [TTj) , [2Tj) . 

For the square lattice, the threshold fraction of the random bond percolation 
problem is obtained as 1/2 P^J By using Eqs. (|3-3p and (|3-4p . we obtain 

TD = hT3 ~ L820478 ' p o = l = °- 75 ' 

when J/ks = 1- This is the result for the square lattice. When J/ks = 1, the 
ferromagnetic transition temperature for the pure system, Tc, is 2/ln(l + \/2) (~ 
2. 269) 20 1 . T D does not agree with Tc in this case. The previous numerical results 
are T D ~ l.sP for p = 1/2, T D ~ 1.7CpD for p = 1/2, T FK ~ 1.81 and 1.82^ for 
p = 0.7, and Tfk ~ 1.83^ for p = 0.8. It seems that the obtained result is valid. 

Fig. [1] shows a schematic phase diagram for the ±J model, p is the probability 
that the interaction is ferromagnetic, and 1 — p is the probability that the interaction 
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Fig. 1. A schematic phase diagram for the ± J model, p is the probability that the interaction is 
ferromagnetic, and 1 — p is the probability that the interaction is antiferromagnetic. T is the 
temperature. The paramagnetic phase ('Para'), the ferromagnetic phase ('Ferro'), the Griffiths 
phase ('Griffiths Phase') and the spin glass phase ('Spin Glass') are depicted. The Nishimori 
line (the dashed line), the Griffiths temperature (the dotted line) and the Td line (the short 
dashed line) are also depicted. Tc is the Curie temperature for the ferromagnetic model. The 
point 'A' is the dynamical transition point for the distance between two spin configurations on 
the Nishimori line. The point 'B' is the multicritical point. 

is antiferromagnetic. T is the temperature. The paramagnetic phase ('Para'), the 
ferromagnetic phase ('Ferro'), the Griffiths phase ('Griffiths Phase') and the spin 
glass phase ('Spin Glass') are depicted. In this article, we do not mention the exis- 
tence of a mixed phase between the ferromagnetic phase and the spin glass phase. 
The Nishimori line (the dashed line), the Griffiths temperature (the dotted line) and 
the Td line (the short dashed line) are also depicted. Tc is the Curie temperature for 
the ferromagnetic model. The value of the Griffiths temperature corresponds to the 
value of Tq- The Td line represents the dynamical transition temperature for the dis- 
tance between two spin configurations. It has been pointed out that the value of Td 
numerically agrees very well with the value of Tc when p = l.M The phase between 
the Griffiths temperature and the Td line can also be called the Griffiths phase, but, 
in this phase, the behavior of the distance is the same as that in the paramagnetic 
phase. The point 'A' is the dynamical transition point for the distance between two 
spin configurations on the Nishimori line. The point 'B' is the multicritical point. In 
a certain lattice case, the dynamical transition temperature Tp does not agree with 
Tc- On the other hand, in a certain lattice case, the dynamical transition tempera- 
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ture Tjj agrees with Tq. For example, To for the square lattice does not agree with 
Tq for the same lattice as mentioned above. On the other hand, Tp for the Bethe 
lattice agrees with Tq for the same lattice as mentioned above. For the multicritical 
point, it seems that the temperature at the multicritical point generally does not 
agree with Tp. For example, when J/ks = 1, the temperature at the multicritical 
point for the square lattice is roughly equal to 0.957, 12 ^ while To for the square 
lattice is 2/ In 3 ( ~ 1.82) from the present result. In addition, the temperature at 
the multicritical point for the Bethe lattice is 2J/ks m(\/5 — 1 + 1/y/z — 1 — 1) from 
the result in Ref. 23), while Tjj for the Bethe lattice is 2J/ks h\.(z/z — 2) from the 
present result. 

For the simple cubic lattice, the threshold fraction of the random bond percola- 
tion problem is numerically estimated as Pc ~ 0. 248813.® By using Eqs. (|3-3p and 
(|3-4p . we obtain 

T D ~ 3.9347 , p D ~ 0.62441 , 

when J/ks = 1- This is the result for the simple cubic lattice. The previous 
numerical results are T D ~ 4.]P for p = 1/2, T D ~ 4.CP for p = 1/2, To ~ 3.92^3 
for p = 1/2, and T FK ~ 3.9# for p = 1/2. It seems that the obtained result is valid. 

For the 4-dimensional hypercubic lattice, the threshold fraction of the random 
bond percolation problem is numerically estimated as Pc ~ 0.16013P3 By using 
Eqs. (|3-3p and (|3-4p . we obtain 

Td ~ 6.191 , p D ~ 0.5801, 

when J/ks = 1- This is the result for the 4-dimensional hypercubic lattice. The 
previous numerical results are T D ~ £3 for p = 1/2, and T D ~ 6.03 211 for p = 1/2. 
The obtained result does not contradict with the previous results. 

For the triangular lattice, the threshold fraction of the random bond percolation 
problem is obtained as 2 sin(7r/18)P^ By using Eqs. (|3-3|) and (|3-4|) . we obtain 

2 l + 2sin(vr/18) 

T D = , ,„ . , /10 , 2.759641 , p n = ^ — '- ~ 0.6736482 , 

, r l+2sm(7r/18 I ' F o 2 

m Ll-2sin(7r/18)J 

when J/ks = 1. This is the result for the triangular lattice. The previous numerical 
results are T FK ~ 2.74 and 2.73^ for p = 0.6, and T FK ~ 2.75 and 2.76^ for 
p = 0.7. It seems that the obtained result is valid. 

From these results, it seems reasonable to consider that the present theory is 
roughly good, although more precise numerical confirmations are needed in order to 
decide the correctness of the present theory. 



§5. Concluding Remarks 



We theoretically showed the conjectured exact locations of the dynamical tran- 
sition points for the ±J Ising spin glass model based on a percolation theory and 
conjectures. 

The present theory may be applied directly to the Gaussian Ising spin glass 
model and the Potts gauge glass model. Concretely, the solution of [{Pfk)t]r on 
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the Nishimori line for the Gaussian Ising spin glass model is obtained in Ref. fl"5|) . 
and the solution of [{Pfk)t]r on the Nishimori line for the Potts gauge glass model 
is obtained in Ref. I27|) . By using these solutions instead of the solution for the ± J 
Ising spin glass model, the conjectured exact equations in the Gaussian Ising spin 
glass model and the Potts gauge glass model are obtained. 

We studied the transition for the time evolution of the distance between two 
spin configurations on the Nishimori line. This study is different from the study of 
the aging phenomena on the Nishimori line as in Ref. [5]) . 
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